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Abstract 

For three anyons confined in a harmonic oscillator, only the class of states 

that interpolate nonlinearly with the statistical parameter contributes to the 

third virial coefficient of a free anyon gas. Rather than evaluating the full 

three-body partition function as was done in an earlier publication (Phys. 

Rev. A46, 4693 (1992) ), here only the nonlinear contribution is calculated, 

thus avoiding delicate cancellations between the irrelevant linear part and the 

two-body partition function. Our numerical results are consistent with the 

simple analytical form suggested recently by Myrheim and Olaussen. 
PACS number(s): 05.30.-d, 03.65. Ge, 05.70.Ce, 74.65.-Fn 



1 



Typeset using REVT^ 



In this note we report an improved numerical calculation of the third virial coefficient of 
a free anyon gas using the method of the hyperspherical coordinates of an earlier publication 
||1|]. The virial expansion of the two-dimensional anyon gas may be written as 

^ = 1 + A2{a){pX^) + A,{a){p\^f + (1) 

where P, p and (3 are the pressure density and the inverse temperature in units of the 
Boltzmann constant, and A = /i(2'7r/3/m)^/^ is the thermal wavelength. The expansion 
parameter is the dimensionless quantity (pA^), and ^42,^3, are the second and third virial 
coefficients of the g function of the dimensionless statistical parameter a. Note that 

the coefficients ^2,^3 etc. are temperature independent, since there is no length scale in 
the statistical interaction. To evaluate the virial coefficients it is convenient to confine the 
anyons in a two-dimensional harmonic oscillator of frequency a;, and take the limit of the 
appropriate combinations of the partition functions as {huj(3) — > 0. Denote the N-body 
partition function in the harmonic oscillator by Z^. The centre-of-mass part Zi may be 
separated out, and the relative part is given by Z^ = Z^/Zi. Then the second and the 
third virial coefficients are given by {x = hujf3) 

A2{a) = lim A2{a,x) 

= limx-2[l-2Z2/Zi], (2) 
A-i{a) = lim A3{a,x) 

x^O 

= hm 2x-^[l + 8(Z2/Zi)2 - 5(Z2/Zi) - 3{Zs/Z^,)]. (3) 

The calculations are performed in the fermionic basis, i.e. for a = 0, the particles are 
noninteracting fermions while for a = 1 they behave like noninteracting bosons. In Eqs.(2) 
and (3), Zi and Z2 are analytically known {t = exp{—x)): 

= - ^)^ Z2 = + - t'f, (4) 

while Z3 in Ref. was calculated from the numerically obtained spectrum of the three- 
anyon system. Three particles in two dimensions, after elimination of the centre-of-mass, 
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have four independent degrees of freedom. The problem may be regarded as that of one 
body in 4-dimensions. Moreover, it is convenient to use the hyperspherical coordinates 0] 
of a 4-dimensional sphere with radius R and three independent angles 6, (j) and ip. This is 
because in these coordinates, the fermionic (or bosonic) basis of states may be constructed 
by specifying their quantum numbers using a simple set of rules . The three-anyon 
Hamiltonian may be separated in a radial and an angular part, and the eigenvalue spectrum 
obtained by diagonalizing the angular part of the Hamiltonian in the basis of the 
hyperspherical harmonics YN^u^x{(^,4','ip)- In addition, there is another quantum number 
n' associated with the radial excitations. The eigenstates of the harmonic oscillator are 
completely specified by the set of quantum numbers (A^, u, A, n') with eigenenergies given by 
EN,n' = {N + 2n' + 2)huj. The form of the Hamiltonian and other mathematical details will 
not be repeated here, since these are given in Ref. 

The calculation of in Ref. involved a delicate cancellation of terms in the square 
brackets on the right-hand side of Eq.(3). Since ^3 is known to be finite there should 

be no terms of order x and within the bracket, and all terms of order should cancel 
exactly as —* 0, if the numerical calculations were exact. In practice, due to the truncation 
of the basis and numerical inaccuracies, these are not satisfied. The numerical calculations 
in Ref. [||] were reliable for > 0.8, and careful interpolation was necessary to reach the 
limit of x^ = 0. Part of the error resulting from this incomplete cancellation takes place 
between the terms involving Z2 and only that portion of Z3 which originates from the linearly 
interpolating eigenvalues of the three-body problem. The analytical expression of this part 
of Z3 is known [^,|^ . We may decompose Z3 into two parts: 

z, = # + Z3^^ (5) 

where 

=^ ~ (t- 1)2(1 + t + t2)2' W 

is the portion from the linearly interpolating states, and Z^^ is the nonlinear part. Although 
Z^^ diverges as x~^ for x — 0, the quantity 



AZ^^{a) = [Z^\a)-Z^'m, (7) 

is finite. Using Eqs. (3)-(6), it is straight forward to sliow tliat 

Av43(a) = A^{a) - ^(O) = -6 lim ^Z^^{a). (8) 

It follows from the above equation that it is only necessary to compute ^Z^^[q.) numer- 
ically. This necessitates identifying and excluding the linearly interpolating states (after 
diagonalizing the three-body Hamiltonian) from the calculation of the partition function. 
These linearly interpolating states change by 3 units of huj as the statistical parameter a 
varies from to 1, and may be identified by the set of rules given in the Appendix. 

The diagonalization of the angular part of the Hamiltonian (as specified in Ref. was 
performed with a basis truncated at 

Nmax= 130 for < |A| < 20 

= 100 for 21 < |A| < 72. (9) 

The quantum number A determines the angular momentum of the three-body system. Note 
that this basis for diagonalization is considerably bigger than used in Ref. IQ. After the 
diagonalization, only the nonlinearly interpolating states were used to calculate /S.Z^^[q.) 
as a function of for various values of a. The nodeless {n' = 0) states were included 
in the evaluation of X^i exp(/3i?j), and the infinite tower of nodal excitations were taken 
into account by the multiplicative overall factor of [1 — exp{—2x)]^^. The numerical results 
changed insignificantly even if the sum over all the eigenstates, including the tower of nodal 
excitations were cut off at 90huj. The quantity [—6AZ^^{a)] is computed as a function 
of for a = 0.05 to 0.5 in steps of 0.05. In principle, its value at a;^ = should yield 
AAs{a). In practice, there is still a divergence due to numerical inaccuracies, since the 
cancellation of singularities between Z^^{a) and Z^^{0) is not perfect as x"^ —>■ 0, since 
Z^^{a) is calculated less accurately than Z^^{0). The spurious singularities in AZ^'^^a) 
have to be isolated to extract Ay43(Q;). To this end, we define AA^^a, x^) = —6AZ.^'"{a) for 
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any x, and fit ^ this function by the following form containing six a-dependent parameters 
in the range 0.05 < < 3.4: 



a 

The computed AZ^^{a) is sensitive to the basis size only for < 0.2. 

In the above fit, the coefficient c = AAs{a), and is found to be numerically stable even 
when higher order terms biQX^^ + bi2X^'^ + 6i4X^^ are added in the exponent of Eq. (10). The 
parameter n in Eq.(lO) is close to 2 upto a = 0.20, but decreased monotonically for larger 
a, with n = 1.56 for a = 0.50. The convergence of the extracted AA^{a) values are also 
tested by performing the same calculation of AZ^^{a), but including only eigenvalues upto 
Njnax = 60. The Ay43(Q;) values so obtained are displayed in Table 1 for various values of a 
in the range 0.05 to 0.40. For a > 0.40, the extracted AA^la) is not too stable, and is too 
inaccurate to be listed. This is because the eigenvalues of the Hamiltonian matrix for the 
larger a values tend to be more and more inaccurate. 

Numerically, of course, only A3 values upto a = 0.5 would suffice, because of the relation 
Asia) = ^3(1 — a) derived by Sen ^ . Recently, Myrheim and Olaussen P have calculated 
A3 (a) numerically using a path integral representation of the three body partition function. 
Their method of calculation is completely different, and uses the winding number formalism 
rather than the statistical interaction. From their extensive numerical work, they suggest 
that perhaps Ay43(Q;) obeys the remarkably simple relation 



AA3{a, x'^) = — + c expffeax^ + hx'^ + b^x^ + bsx^]. (10) 
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AA-sia) = — (11) 

For small a, this is in agreement with the perturbative analytical result of a^/12 of de 
Veigy and Ouvry 0. In Fig. 1, the computed values AA^i^a) of the present calculation are 
compared with the simple form given by Eq. (11). Our calculated Ay43(a;) values, confirm 
Eq. (11) with good accuracy upto a = 0.25, beyond which our numerical method starts 
to lose accuracy. This independent check of the Myrheim-Olaussen calculation is the main 
content of this note. Fig. 1 also shows that the inaccuracies for larger a build up rapidly in 
the method of Ref. [jl|, where the nonlinear part Znl was not isolated. 
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APPENDIX 

The exact solutions for the linear states with slopes of ±3huj have been written down by 
Wu [0 in Jacobi coordinates. Using the details given in Ref. f^, we can easily write down 



these solutions in hyperspherical coordinates and identify the quantum numbers A^, z/, A, n' 
and the degeneracy of these states. We find that for A^ = {6p — 3), {6p — 1), or (6p + 1) 
(with p = 1,2,3...), there are degenerate states with n' = 0, (and hence total energy 
E = N + 2), while for A^ = 6p, {6p + 2), or {6p + 4) (with p = 1, 2, 3, ...) there are 2p{p + 1) 
degenerate states with n' = 0. For example, for A^ = {6p — 3), one finds that the states with 
A = ±(6p — 3) are both p— fold degenerate with the corresponding u quantum number being 
u = (6p — 3), {6p — 9), ...3 (note that z/, A^ > 0; u, |A| < A^; u, N, |A| have the same parity, and 
u is always a multiple of 3). For the same A^, the four states with A = ±(6p — 5), ±{6p — 7) 
are all (p — 1)— fold degenerate = {6p — 3), (6p — 9), ...9). The next four states with 
A = ±{6p — 9), ±{6p — 11) are all {p — 2)— fold degenerate (z/ = {6p — 3), {6p — 9), ...15), 
and so on. Finally the four states with A = ±(2p + 7), ±(2p + 5) are all non-degenerate. 
Exactly the same rules also apply when A^ = {Qp — 1) or {6p + 1). On the other hand, when 
A^ is even, for A^ = 6p, we find that the four states with A = ±6p, ±(6p — 2), are all p— fold 
degenerate with u = Qp, (6p — 6), ...3. The next four states with A = ±(6p — 4), ±(6p — 6) are 
all (p— 1)— fold degenerate and so on. Finally the four states with A = ±(2p + 2), ±(2p + 4) 
are all non-degenerate. Similar rules also apply when A^ = (6p + 2) or (6p + 4). In this 
way we are able to identify all the linear states in the spectra and exclude them from the 
evaluation of the partition function. 
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FIGURES 

FIG. 1. Plot of AA3{a) as a function of the statistical parameter a in the fermionic basis, [a] 
Present calculation (top row of Table 1), with the six-parameter fit of Eq. (10). [b] The Variations 
in the extracted A^3(a), with a nine-parameter fit to AAs{a,x'^). This is done by including 
the terms 6102;^° -|- bi2X^'^ + hux^^ in the exponent of Eq. (10). [c] The plot of Eq. (11), (the 
Myrheim-Olausscn conjecture for /S.A'^{a)). [d] The calculated lS.A^{a) values by the method given 
in Ref. [1], but in the enlarged basis, Eq. (9). 
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TABLES 

TABLE I. The calculated quantity AA3(a) = [A3 (a) — A3(0)] is tabulated as a function of a. 
The first row of A A3 (a) is obtained by including all the cigcnencrgies in the space defined by Eq. 
(9). The six-parameter fit, as given by Eq. (10), is used. The second row is computed by including 
only eigenvalues upto Ecut = 60^a;. 



a 


0.05 


0.10 


0.15 


0.20 


0.25 


0.30 


0.35 


0.40 


AA3(a) X 10"^ 


2.08 


8.22 


17.3 


29.0 


41.8 


53.8 


62.7 


65.7 


AA3(a) X 10'' 
Ecut = QOhio 


2.13 


8.34 


17.6 


29.4 


42.4 


54.7 


63.9 


67.4 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9310106vl 



